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ABSTRACT 


This Report contains the results of a one year study into the 
effects of a pulse shape on the transient response of a cylindrical 
shell. Uniaxial, membrane, and bending theories for isotropic 
shells were used in this study. In addition to the results of the 
above analytical study, the preliminary res ul ts of an experimental 
study into the generation and measurement of shear waves in a 
cylindrical shell are included. 
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NOMENCLATURE 


K 2 = shear correction factor 
c 2 = plate velocity = — 

p P (w) 

/■» 

= shear velocity = K - 

p 

t Q = pulse duration 

t C 

t = dimensionless time = 
t q - dimensionless pulse duration 
X = equivalent pulse length = Cpt 
X = dimensionless pulse length = 

Other symbols are defined in text 



I. INTRODUCTION 


Included in this Report are the results of a one year study performed 
for the Structures Division of NASA - Langley Research Center. The major 
tasks of this study were: 

1. An analytical parametric study to determine the effect of pulse 
shape (e.g. magnitude, shape, duration, rise time) on the transient 
response of cylindrical shells subjected to longitudinaT impacts, 

2. An analytical parametric study to determine the effect of pulse 
shape on the transient response of cylindrical shells having geometri 
cal discontinuities (e.g. discontinuity in thickness) subjected to 
longitudinal impacts, 

3. An evaluation of the importance of secondary shell theory terms 
(e.g. transverse shear deformation, radial and rotary inertia, and 
bending) in predicting transient responses of cylindrical shells 
subjected to axial impacts. 

In addition to the above tasks, an experimental program was initiated 
to develop a technique for generating shear waves in a cylindrical shell. 
The initial results of this program, including the measured shear wave 
velocities, are also included in this Report. 



II. SOLUTION OF THE MEMBRANE SHELL EQUATIONS 
BY THE METHOD OF CHARACTERISTICS 


Since part of our task in this grant was to determine the effect of 
the higher order shell terms (eg rotary inertia, transverse shear strain, 
bending) on the predicted responses of cylindrical shells subjected to 
longitudinal impact, we decided to attempt the solution of the classical, 
membrane equations by the method of characteristics. The first point to 
be mentioned is that the classical membrane equations do not constitute a 
system of completely hyperbolic partial differential equations. The method 
of characteristics is only applicable to hyperbolic systems of equations. 
Obviously the application of the method of characteristics to a non-hyper- 
bolic system of equations is suspected. Before discussing the numerical 
approach used here to achieve this apparent misapplication of a mathematical 
technique, we shall describe the classical membrane equations. 

The equations of motion of a cylindrical membrane under axisymmetric 
conditions are (see Appendix A for equations of all shell theories used) 

3N X a 2 u 

3X " ft 7 

n 9 a 2 ,, , m-D 

IT ‘ hp 


and the strain-displacement and stress-strain relations are 
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where N , e , and u are the stress, strain, and displacement in the 

X XX 

axial direction; N. and e aa the stress and strain in the circumferential 
directions; w the displacement in radial direction; h and R the shell 
thickness and radius. 

Substituting equations (II-2) into (II-l) yields the classical 

t 

membrane equations in terms of the two displacements. These are 

3 2 u , a 2 u 3 w 

1 _ v_ 

3X Z ~ Cp 2 dt T R 3X 

(II-3) 
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T^v 2 ^R 


3U 

+ V 3^r ] 


where Cp 2 = E/p(l-v 2 ) called the plate velocity. As we see, this system 
of equations is not completely hyperbolic. The first equation can be 
considered to be hyperbolic (the left-hand side is of the form of the simple 
wave equation) while the second equation is parabolic in nature. In order 
to understand our reasoning, soon to be introduced, let me describe another 
"membrane" theory which will incorporate the effect of the transverse shear 
force into the classical membrane formulation. For this case the "membrane" 
equation of motion are (Ref. 1) 


5 

3X 


= ph 


3 2 U 

St 2- 


!!k Vph 3 !^ 

3X ~ R 


3t^ 


(n-4) 


where in addition to the relation of ( I 1-2) we how have 


3 W 


3 W 


Y XZ 3X and ^X K2y XZ ^ G 3X 


(H-5) 


where. G is the shear modulus, Q the transverse shear stress, y the trans- 

x xz 

verse shear strain, and K 2 the shear correction factor. Substitution of 
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equations ( I I -2 ) and ( I 1-5 ) into ( I I -4} yields the system of equations 

3 2 U , 3 2 U 3-W 

- _ I . _ v_ 

3X 7- ~ Cp 2 3t 2 R 3x 

(I 1-6) 



where c g 2 = K 2 G/p is called the shear velocity. Equation ( 1 1 -6) is a com- 
pletely hyperbolic system of partial differential equations; disturbances in 
u will propagate with the plate velocity, Cp, and disturbances in w will 
-propagate with the shear velocity, c $ . We can now see, conceptually, the 
difference in physical interpretation between equations ( I I -3) and ( I I -6) . If 
in the second of equations (1-6) we multiply through by c s 2 (or K 2 ) and then 
let c s 2 go to zero we see that eq. ( I I -6 ) reduces to eq. ( I I -3) . So, our 

i 

numerical procedure used here will be to solve eq. (II-6), but, we will 

require c g 2 to be extremely small, but , not zero. In other words, we will 

actually be solving the systems of equations ( I 1-6) , but, due to the smallness 

of c s 2 we are physically approximating equations ( I I -3) . The question now 

arises, the application of the method of characteristics to equations (I 1-3) 

has been a standard technique so why this procedure? The answer is simply this. 

When people do apply the method of characteristics to equations (I I -3) in order 

to analyze a cylindrical shell impact problem they always consider a semi - 

infinite medium . The reason for this is simple, this standard technique does 

not permit the incorporation of a boundary condition in w (remember a term 
3W 

appears in the governing equations). So, at the impacted end a boundary 
condition in u is prescribed and the other end is assumed to be at infinity, 
thus, the question of imposing one boundary condition in w is avoided. With 
our procedure outlined here we are able to mathematically incorporate all the 
boundary conditions properly. The only point to be demonstrated is whether 
oun solution of equations (1 1-6) with c s (K ) very small yields, for practical 
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purposes, the solution to the classical membrane equations. 

We ran some test cases to numerically determine how small c g 2 must be. 
Our procedure was to use MCOIT-21 (Ref. 2) to solve equations (II -6 ) 
subjected to impact boundary conditions at x=0, or, 

3U ’ 

at- (°’t) - v 

3 W 

G K 2 (0,t) =0 

We used the shear correction factor, K 2 , as our parameter to vary the shear 
velocity, c s (remember c g 2 = K 2 G/p); the value of K 2 is approximately 0.87 
for most shell problems involving shear waves. We calculated the transient 
responses of a typical cylindrical shell subjected to axial impacts with 
different values of K 2 . We then compared the strains at a particular point 
as predicted by each run. When the strain predictions did not significantly 
change with a further reduction in K 2 we chose that value of K 2 for our 
later studies involving the classical membrane theory. As an indication of 
our K 2 study, Table (II -1) shows the comparison of longitudinal strain at 
a distance 60 thicknesses from the impacted end for various values of K 2 . 

TABLE (II -1) 


K 


0.87 

.75649 

0.87 x 10" 4 

.77251 

0.87 x 10" 6 

.77318 

0.87 x 10“ 8 

.77391 

0.87 x 10' 10 

.77225 

(double precision) 
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We can see from Table (II -1) that, numerically, once K has reached the 

-4 - 

value of 0.87 x 10 there is no significant change in the strain predic- 

2 

tions with a further reduction in K . Except for the last entry all other 
calculations were performed in single precision. 

In order to determine the exactness of our technique for solving the 
classical membrane equations we compared our solution to the Laplace trans- 
form solution of the classical membrane equations as published by Berkowitz 
(Ref. 3). ‘ In his paper, the author analyzes the response of a semi- 
infinite elastic cylindrical shell subjected to a longitudinal impact. He 
solves a system of equations identical to (II -3), subjected to impact 
conditions at x=0, by applying asymptotic expansion techniques to the Laplace 
transform inversion integrals. He calculates the longitudinal stress at 
a dimensionless time as it varies with dimensionless axial distance. Figure 

2 _o 

1 shows ... his results as compared to our results for K = 0.87 x 10 . 

2 o 

The fact that our results (for K = 0.87 x 10" ) and h is do not agree identi- 
cally is to be expected since both are approximate solutions to the same 
governing equations. 
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III. PULSE SYNTHESIS FOR- DETERMINING TRANSIENT. RESPONSE OF SHELL DUE TO , 
PULSES OF DIFFERENT SHAPES 

A major portion of this research project was to determine the effect 
of pulse shape (e.g. rise time, pulse duration, shape) on the transient 
response of cylindrical shells subjected to longitudinal impact loading. 

One approach for achieving this phase of the research is to utilize our 
computer codes (MCDIT-21 and MCDU-26) and make numerous runs where we vary 
the shapes, rise times, and pulse durations. The result of this type of 
approach would be numerous parametric plots. Another approach would be to 
determine a set of "building blocks" which could be used to approximate the 
different pulse shapes and then we need only evaluate and understand the 
response of the shell to these "building block" functions. Once we under- 
stand this then we need only superimpose the resulting responses linearly 
and we can understand the response to the original pulse shapes of interest. 
This letter concept is not new; people have discussed the use of rectangul ar 
functions as "building blocks" for such a purpose. However, we discarded 
this type of function for our analysis due to the large number which would 
be needed to describe the pulse shape functions used in our study. The 
approach we finally chose to achieve this phase of our research was a 
combination of both. In otherwords, we performed some computer runs with 
the exact shapes in order to isolate important parameters of the pulse 
and then we used the second approach to try and understand these initial 
results and predict responses with further variations of the parameter. 

In this Section of the Report we will discuss the pulse shapes used in our 
study, the. "building block" approach, and the response of the cylindrical 
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shell to these "building block" pulses. In the next Section of the 
Report, we will, discuss the exact approach, the use of the "building 
blocks", and the results- of our pulse shape study. 

III. 1 Pulse Shapes 

We decided to limit the number of pulse shapes to those shown in 
Figure 2. The reason we chose these shapes is that we felt that each 
shape or combination of shapes was of practical interest. For example, if 
we understood the response of the shell to shape 2a and 2b we could then 
understand the response to a pulse whose rise time ranged from 0 to t q /2 
(see Figure 3a). This latter shape can be seen to closely represent a 
typical explosive pulse (Figure 3b). Other practical pulse shapes can 
be seen to be composed of the shapes shown in Figure 2. One last point 
should be made here. Each of the shapes shown in Figure 2 yield 
identical impulse values so that when we compared the responses due to 
these shapes there was- no difference in energy input. 

III. 2 "Building Block" Synthesis 

The essence of this principle, for linear differential equations, is 
simply as follows: 
if 

a(Input A) + b(Input B) = (Input C) 

then 

a(Solution for A) + b(Solution for B) = Solution for C 
where a and b are linear coefficients. In otherwords, if we -can represent 
a particular function by a linear combination of approximating functions 
then we can approximate the solution for the original function by super- 
imposing the solutions for the approximating functions. 
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Before describing, in detail our search for these approximating 
functions {"building blocks") we should define a few terms used in the 
discussion 

Zero order discontinuity - discontinuity in magnitude only 

First order discontinuity - discontinuity in first derivative 

only 

Nth order discontinuity - discontinuity in Nth order deriva- 
tive only 

Approximate linear relation is 

F(t) = A 1 V 1 (t) + A 2 V 2 (t) + A 3 V 3 (t) .... 

where 

F(t) - approximated function 

( t) , V 2 (t), V 3 (t) - approximating functions 

Ai , A 2> A, - linear coefficients 

In our search for the approximating functions we had two obvious 
requirements. First, these functions must be able to approximate the five 
shapes of interest (Figure 2) and, second, the deviation between the exact 
shapes and the combination of these approximating functions should be small 
with only a few number of these functions being used. Once the approximat- 
ing functions were determined we used a least squares technique to determine 
the linear coefficients Ap A 2 , etc. A computer program was written to 
facilitate the determination of these coefficients for any shape approximation 
In addition, this program plotted the approximated shape based on the number 
of approximating functions desired. This computer program is contained in 
Appendix B. 

Described below are some of the approximating functions considered 

A. Fourier Series - Single Pulse 
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The initial specifications for the approximating functions was that 
they should have no zero order discontinuity at the origin and they should 
be continuous from 0+ to + » ( should not have a discontinuity of any 
order in the 0+ to + ® region). The above requirements are immediately 
satisfied- .by' the sine series, (sin nut ; n »= 1,2,3 ...) The cosine 
series was not considered because of the zero order discontinuity on the 
origin. Due to the nature of the approximation and due to the fact that 
the sine series is periodic, one single pulse can not be approximated. 

For example; when the square pulse is approximated what in effect has been 
created is a new periodic function consisting of a series of square pulses. 

B. Fourier Series - Compound Pulse 

Observing the results of the previous example we see that an approxi- 
mation of a pulse by a continuous periodic series produces another series. 

To partly alleviate the above problem at new pulse is created, which is 
composed of the original pulse plus a finite zero region. The Fourier 
series approximation produces again a new series, but it is hoped in this 
case, if the zero region length is quite large with respect to the square 
pulse length, and the approximation is reasonable, for all practical 
purposes it will mathematically satisfy the requirements of a pulse. 

This compound pulse however can not be easily approximated. The 
results shown in Figure 4 indicate that even after ten terms the approximation 
was still quite poor. 

C. Fourier Integral 

The pulses can be approximated by continuous periodic functions not 
by a series, but rather by a continuous integral. This is of course useless 
for the study because no discrete functions can be extracted. 
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D. The Envelope Function 

Certain continuous functions can form a pulse of a variety of shapes 

as was shown by the Fourier integral. The Fourier integral though is 

difficult to evaluate. Another function that forms a square pulse is 
-t m 

e As "m M approaches + » function becomes square 
This envelope function can be used as a multiplier for other Fourier series. 
The use of the envelope is to accentuate the part of the periodic series 
which represents the pulse, 

E. Raising Sine to a Power 

For the sake of interest the sine series was raised to a reciprocal 
power and its approximating ability was observed. This widening effect 
was excellent for approximating the square pulse but it failed on the 
triangular pulses. 

F. Creating a Function 

In our study we are limited to only the use of three approximating 
functions. This means that the approximations will be poor especially 
for the type of functions of interest. To alleviate this problem a 
special set of approximating functions are created. These functions can 
be made to form an almost exact fit to the pulses under consideration. 

One effective way of accomplishing this is by replacing some sine functions 
by polynomials. 

G. Removing Unwanted Portion of Fourier Functions 

From part "A" we see that the Fourier series produces periodic functions 
instead of pulses. Instead of using the complete function, we will only 
utilize the portion up to t . For the sine series however, by removing 
the unwanted portion of the periodic function you introduce discontinuities 
of order one and higher. The utilization of this technique to represent our 
pulse shapes is demonstrated in Figures 5 (three terms) and 6 (five terms). 


- 15 - 




FlS Otel- 5 ** TUR.E E APPKova 






















SUMMARY 


1. A continuous set of sine functions cannot be used since the 
approximation creates another periodic function instead of a 
single pulse. 

2. A Fourier integral cannot be used since it is not composed of 
a sum of discrete functions. 

3. The envelope can be used but it is necessary since first order 
discontinuities are not harmful for computer applications. 

4. Raising a sine to a power increases the deviation for some of 
the shapes. 

5. Replacing some of the sine functions by polynomials is only use- 
ful for obtaining accurate approximations to only a limited 
number of pulses. 

6. Using the sine series and removing the portion of the functions 
after t q can be used. The only problem here is due to the first 
and higher order discontinuities which are introduced by removing 
part of the function. For computer applications, however, this 
problem is negligible. 
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Having decided on the last approximation technique we decided to 
compare the results of a three term approximation for a rectangular 
pulse with the results obtained for the exact shape. The problem we 
chose was that of the longitudinal impact of a cylindrical aluminum shell 
(h/R = 0.049). The boundary conditions for this problem were M =Q = 0 

A A 

and 

0(0, t) = 75 in/sec ; 0£t<_ 21 ysec 

N x (0,t) =0 ; 21 ysec <t 

In otherwords, the pulse shape was rectangular. We used the three term 
approximation for the rectangular pulse (see Figure 5). We then solved the 
governing system of equations, eq. Arl , for the bending theory (k- - .87) 
for each of the three approximating sine functions. The longitudinal strain, 
e x , at X=0 and 3 inches for each of these functions is shown separately in 
Figures 7a and b. We then superimposed the response of each function 
(by use of the linear coefficients of Figure 5) and compared the resulting 
response with those of the numerical solution using the exact shape and 
experimental results. The longitudinal and circumferential strain compari- 
sons are shown in Figures 8a and b. We see that the three terms approxima- 

* 

tion solution agrees quite well with the solution for the exact shape. 

III. 3 Transient Response of the Sine Approximating Functions 

Our next step in the "building block" approach was to evaluate the transient 
response, of a cylindrical shell to the sine pulse, the. function chosen as our 
approximating function. Of course,, by .following this approach of determining the 
effect of pulse shape on a shell's response we eliminate a parameter, namely, the 
actual shape of the pulse (including rise time). Our remaining parameter is the 


Both of the solutions were based on the bendina theory, eq. A-l 
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pulse duration. The problem now reduces to determining the effect of the 
duration of the sine pulse on the transient response of the shell. If 
this is understood, then by superposition we should be able to determine 
the effect of an arbitrary pulse's shape and .duration on the shell response. 

Our first set of runs involved a cylindrical shell (h/R =0.1) 
subjected to an axial velocity input of the type 

u(0, t) . = sin rr ~ 0< t< t 

T — — 0 

0 

N x (0, t) = 0 t q <t 

where t q is the pulse duration of the half sine. The remaining boundary 
conditions at the impacted end were Q =M = 0. Transient responses, as 

A A 

predicted by the "bending theory" eq. (A-l ), were obtained for values of 
t q such that the inverse of wavelength, 1/A = l/x o (or in dimensional form A = 
t Q c p I assumed values of .2, .1 , .05, .025, .0167, and .01. In otherwords, the 
range of the pulse duration was such, that the equivalent pulse length varied 
fronj, being 5 times the thickness to TOO times the thickness-. Plots of the 
longitudinal strain at a location of X = 60 are shown in Fig. 9. Comparison 
of these responses demonstrates an interesting point; the sine pulse disperses 
less for the long and the short pulse duration than for the intermediate 
duration. Responses for identical loading conditions were obtained for the 
.classical membrane theory eq. (A-3), the membrane with shear theory eq. (A-2), 
and the uniaxial theory eq. (A-4) . The trends of the longitudinal strain as 
predicted by the first two of these theories were identical to that of the 
"bending" theory although there was a slight deviation in magnitude. Of 
course the uniaxial theory indicated no dispersion across the entire wave- 
length range since it is governed by the simple wave equation. To illustrate 
this pulse length effect on the transient response of the cylindrical shell 
we plotted in Fig. 10. the peak values of the longitudinal strain versus the 
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(or at the location. X ='60. A similar plot of this variation for 
a cylindrical shell with h/R = 0.05 is shown in Fig. 11. This variation 
of strain with pulse length can be described subjectively as follows: 

L. 

when the wavelength is large 0), as compared to the R of the shell, 
we can consider the wave as seeing a bar and it is well known that for 
a wave length much larger than the bar radius there exists no dispersion. 

On the other hand, when the wavelength is smaller than the shell radius 

L 

and of the order of the shell thickness > 0.5) the wave is seeing 
essentially a plate. For this case we can understand the lack of dispersion 
(mathematically) by realizing that the inplane equation of motion 
(uncoupled from the bending and transverse equations) for a plate is 
governed by a simple wave equation. The region between these two limits 
involves coupling between the membrane, radial, and bending modes, thus 
the dispersion. 

The essence of this discussion, as it pertains to understanding the 
effect of pulse shape and duration on the transient response, is as follows. 
First, I believe we have demonstrated the importance of pulse duration 
(equivalent wavelength) on the transient response of a shell. As we have 
seen, it is the most important parameter to be considered. Second, by 
using the "building block" approach we can understand the effect the pulse 
shape will have on the response by simply determining the importance of 
each of the approximating functions (simr.t, sin 2tt£, sin 3Tit, etc), for a 


particular shape via the least square coefficients. Once we know which of 
these functions is most important we can use curves such as those of Figs. 

10 and 11 to determine the magnitude and dispersive nature of each and thus 
be able to predict the response of actual shape through superposition. This 
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latter point will be discussed in more detail in Section IV where we will 
compare transient responses as obtained by the "building block" approach 
and the exact shape approach. 

Before concluding this Section, a few points should be mentioned. 

First, work should be continued along the line of effort just described 
as I believe this will lead to the understanding of the relationship 
between dispersion analysis (harmonic) and actual wave propagation. Second, 
this understanding of the effect of pulse duration on dispersion for 
single-type pulses has importance in ultrasonic NDT work (see for example 



IV. EFFECT OF PULSE SHAPE ON TRANSIENT RESPONSE OF CYLINDRICAL SHELLS 

In this Section of the Report, we demonstrate the importance of the 
pulse shape (especially the pulse duration) on the transient response of 
a cylindrical shell. In achieving this goal, we will compute the response 
of the shell due to the exact shapes used for the loading functions. How- 
ever, in order to demonstrate the usefulness and importance of our "building 
block" concept, we will choose certain of these shape responses and compare 
them with the predictions obtained by this "building block" approach. By 
doing this, the importance of certain of the pulse shape parameters can be 
isolated and described. In deciding upon which shapes, specimen properties, 
and pulse durations to use in these final computations, we relied heavily 
upon the knowledge gained from the work previously described (e.g. the e vs 
h/X curves of Figs. 10 and 11). This was necessitated by the fact that we 
could become involved in an endless parametric study with the major results 
being lost in the massive data. For the sake of clarity and completeness in 
this Section, certain details previously described in this Report will be 
briefly repeated. 

The two cylindrical shell specimens chosen for our analysis are shown 
in Fig. 12. Each of these specimens is seen to have a discontinuity in 
thickness. The reason for this is that a portion of the proposed study per- 
tains to the effect of geometrical discontinuities (the results are discussed 
in the next Section), so, by analyzing the specimens, as shown, we then have 
data necessary for the uniform shell and the geometrical discontinuous shell. 
Notice that for each specimen L-j/h = 200, or, the length of the shell to be 
considered in this Section is 200 times the thickness. The pulse shapes chosen 
for our study are the same as those shown in Fig. 2. The boundary conditions 
used were 
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J u (0,t) = 1 x shape function; 0 <_ x <_ t q 

n (o,t) = o ; < t 

and 

M x C0,t) = Q x (0,t) = 0 

and the initial conditions were zero. Note that these boundary conditions 
were all applied at the left end of the specimens- as shown in Fig. 12. 

Finally, the equations used for the analysis were those of the bending theory 
(eqs. A-l) since we have seen (Figs. 10 and 11) that, except for the uniaxial 
theory, there is no sizeable discrepancy between the predictions of the 
bending theory with those of the modified membrane theory or the classical 
membrane theory. We should remember two points here. First, the uniaxial 
theory will yield a transient response in which there is no change in shape 
or magnitude of the pulse as it propagates along the shell and second, the 
fact that there is only a slight difference between the predictions of the 
other theories applies because we are only considering longitudinal loadings; 
this is not true for radial loadings. 

The first set of computations involved specimen #1 for which h-j/R-j = 0.05. 
Observation of Fig. 11 shows us that for the pulse shapes having a value of 
1/A = h/x = 0.05 we should not expect much shape change or dispersion as the 
wave traverses the shell. In order to reach this conclusion we are relying 
on our sine "building block" concept when using Fig. 11. Recall that each 
of the exact shapes can be approximated by a series of sine pulses, the first 
sine term having a pulse duration equal to the exact shape pulse duration. 

The second and higher terms of the sine series will have pulse durations 
(or A ' s ) which are shorter, thus, will have even less dispersion because we 
are moving to the right on Fig. 10 (e.g. the third term will have a value of 
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h/x of 0.15). On the other hand, by observing Fig. 11 we would expect 
pulse shapes having a pulse duration equivalent to h/x = 0.01 to suffer 
much dispersion. With these predictions established from Fig. 11, we 
then ran two sets of computer calculations; the first, shown in Fig. 13, 
are the results for the shapes having a pulse duration equivalent to 
h/x = 0.05 and the second, shown in Fig. 14, the results for h/x = 0.01. 

In Figs. 13 and 14 we have plotted the longitudinal strain at x = 50, 100, 

150 and 300; remember that the location x = 300 is beyond the geometrical 
discontinuity, so, the strain history at this location will not be discussed 
until the next section as it represents the transmitted wave through the 
discontinuity. Observation of Fig. 13 and 14 confirms our predictions 
from Fig. 11. For example, in the case of h^/x = 0.05 we note very little 
change in shape, pulse duration, or magnitude as the wave propagates down 
the shell; the peak magnitude of the isoceles triangle pulse at x - 150 
ts 0.91. On the other hand, we see for the case of h-j/x = 0.01 the pulse 
changing shape, increase of pulse duration, and more attenuation of peak 
magnitudes; the peak magnitude of the isosceles triangle pulse at x = 150 
is now 0.75. ■ 

The second set of computations involved specimen #2 for which h-j/R-j = 
0.10. Observation of Fig. 10 shows us that more dispersion is possible 
for this h/R ratio than in the previous case. For this specimen, then, 
we performed three sets of computations; these were for the values of h-j/x 
of 0.1, 0.025, and 0.01. We would predict, based on Fig. 10, that the pulse 
duration causing the least dispersion of the pulse would be for those pulse 
shapes having h^/x = 0.1, while the most dispersive would be for h-|/x = 0.025; 
the pulse shapes having h^/x = 0.01 should lie between these two. Figs. 15, 
16, and 17 are the results of the computer calculations for values of h-j/x 
of 0.01, 0.025, and 0.1, respectively. Observation of Figs. 15, 16 and 17 
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indeed con firms our predictions. The pulse shapes with h^/x = 0.1 did 

not change shape drastically, maintained essentially the same pulse duration 

and had a peak value for the isosceles shape of 0.8 at x = 150. The pulse 

shapes with h-j/x = 0.01 did change shape some what, increased in pulse 

duration, and had a peak value for the isosceles shape of 0.85 at x = 150. 

* 

Finally, the results for the h-j/x = 0.025 show a much more drastic shape 
change and increase in pulse duration than either of the previous two cases; 
the peak value for the isosceles triangle has diminished to a value of 0.65. 
Note, also, from these Figures that the longer pulse duration pulses lose 
their identity in shape and begin to resemble one another in shape as the 
wave propagates along the shell. This indicates that for long duration 
pulses the present practice of not worrying about the shape of the initial 
pulse in detail is most likely proper for practical applications. 

The results of this study indicate that the pulse shape parameter 
which effects the transient response of a cylindrical shell is the parameter 
h/x which is directly related to the pulse duration, r , through the relation 
X = t Q Cp. We have shown that the degree of dispersion of a pulse is related 
to this ratio. A useful tool in predicting whether a particular pulse will 
be dispersive is the type of curves shown in Figs. 10 and 11. In order to 
demonstrate the usefulness of the "building block" concept to actual ly 
predict the transient response of the shell subjected to a particular shaped 
pulse, we will compute the transient response of specimen #2 due to the 
isosceles shaped pulse by this technique. We will use the two term approxima 
tion of the isosceles shape as shown in Fig. 5; in otherwords 

0.81099 sin tt l- - .09048 sin 3n !• 

T T 
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We will perform the calculations for this isosceles shaped pulse for the 
cases of h-j/x = 0.025 and 0.01. The first step is to compute the response 
of the shell to each of the sine functions. Next, we superimpose the 
results by use of the previous equation and we have the resulting approxi- 
mations for the transient response of an isosceles shapes pulse. Figs. 18 
and 19 contain the results of the "building block" response predictions for 
h|/x = 0.025 and h-|/x = 0.01, respectively. These Figures each contain the 
response to the sine functions, the superimposed response, and the response 
due to the actual shape. Comparison of the response predictions as obtained 
by the exact shape and the approximation, show that the "building block" 
concept can be used not only for subjective predictions, but, also practical 
predictions. Once we understand the way in which individual sine pulses 
effect the transient response of a shell we can superimpose and understand 
the response due to an arbitrary shaped pulse. 

In conclusion, we must say that an experimental program to verify these 
analytical results is warranted; the reason being that we are solving a 
system of approximate (although inclusive) governing differential equations 
for the shell response. An experimental verification of the curves shown in 
Figs. 10 and 11 would be sufficient to verify our results. 
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V. EFFECT OF PULSE SHAPE ON TRANSIENT RESPONSE OF CYLINDRICAL SHELLS 


HAVING GEOMETRICAL DISCONTINUITIES 

The purpose of this phase of the research program is to determine the 
effect of pulse shape on the transient response of a cylindrical shell with 
a geometrical discontinuity subjected to a longitudinal impact. In this 
Section we will be referring to Figs. 13 through 19, again, as was mentioned 
previously. In particular, we will be interested in the strain responses at 
x = 150 and 300 in these Figures. The strain histories at these two locations 
can be considered as the incident and transmitted pulses, respectively. 
Observation of Figs. 13 and 14 for specimen #1 and Figs. 15, 16, and 17 for 
specimen #2 establishes certain trends: 

1. For the shorter pulse duration the actual shape is important 
in establishing the peak transmitted pulse magnitude. For 
example, observation of Fig. 17 demonstrates that the maximum 
magnitude is established by initial rise time. 

2. For the longer pulse duration the shape does not play as 
important a role as for the shorter duration pulses. For 
example, the peak magnitudes of the transmitted pulses, as shown 
in Fig. 15, are not dependent upon the rise time or shape 

t 

parameters. As a matter of fact, we see in this figure that 
the peak magnitudes occur well behind the wave front and that 
the deviation between these magnitudes is not as large as that 
which occurs for the shorter duration pulse. Note, also, in 
Fig. 15 that the pulses have essentially lost their identity 
and in fact are quite similar behind the wave front. 





3. No discernible "rule of thumb" formula seems to be appro- 
priate for the prediction of transmitted to incident strain 
or stress ratios. For the shorter pulses having a short 
rise time the simple uniaxial equation (Ref. 5 } 

t £ x] transmitted _ ^1 
[el incident " A 2 +A l 

would yield reasonable results, but, for the pulses with 
either long rise times or long pulse duration it appears 
that a complete computer analysis is necessary. 
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VI. EXPERIMENTAL MEASUREMENT OF SHEAR WAVE VELOCITY 


IN A CYLINDRICAL SHELL UNDER RADIAL IMPACT 

Another phase of research initiated under this grant was that of 
developing an experimental technique for generating and measuring shear 
pulses in a cylindrical shell. The motivation for this study was to 
develop an experimental capability which when combined with our analytical 
capability would equip us with the necessary tools to better understand 
the role of shear waves in the transient response of shells. Following 
is a brief description of this phase of our research effort. 

Experimental studies of propagation of shear wave in cylindrical shells 
are not found in the literature. Some investigators [6,7] studied shear 
wave velocity utilizing the ultrasonic techniques. Steveninck [8] developed 
apparatus for simultaneous determination of longitudinal and shear wave in 
porous media under pressure. He demonstrated the problem of .separating the 
longitudinal and shear wave. The oscilloscope traces obtained in [8] 
demonstrated the difficulty in obtaining a clear strain pulse which propagates 
with the shear velocity. Ref. [?] shows some typical oscilloscope traces 
in the study of radiation from an explosion in a circular disk. 

In this section a technique to measure shear wave velocity in a 
cylindrical shell under radial impact is developed. The concept for 
generating the shear wave was as follows (see Fig. 20). An explosive is 
mounted inside a plug which fits into the end of the shell. The explosive 
is detonated resulting in a pulse propagating radially outward from the 
center. The pulse then induces into the end of the shell a radial loading 
thus generating the shear wave. The radial impact is obtained by an 
explosion of an electrical detonator inserted in plexiglass and in lexan plugs 
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The purpose of using plexiglass or lexan was to slow down the wave initiated 
by the explosion in order that we will obtain a symmetrical wave in the 
aluminum cylindrical shell. A symmetric wave will eliminate the effect of 
bending and thus produce a cleaner strain pulse propagating with the shear 
velocity. 

Strain gages mounted on the outside surface at various axial locations 
of the shell were used to obtain oscilloscope traces; from which the shear 
velocity is determined. However, observation of the oscilloscope traces 
show the presence of a precursor bending wave propagating at the plate velocity 

EXPERIMENTAL PROCEDURE 

A straight cylindrical shell fabricated from 6061 -T6 aluminum was 
radially impacted by an explosion of an electrical detonator embedded in 
plexiglass and lexan plugs. Average properties for the aluminum are 

C 

E = 9.78 x 10 psi = Young's modulus of elasticity 

v = 0.33 = Poisson's ratio 

c = 214,500 in/sec. = Analytical plate velocity = [E/p (1 -v 2 )]^ 2 

r 

c = 113,338 in/sec. -■ = Analytical shear velocity = K(G/p)^ 

V 

( 

The geometrical properties of the. shell are shown in Fig. 21, * 

Two types of electrical detonators were used. The first was the 
Du Pont X-549 R electrical detonator. This type of electrical detonator 
has a slow firing time of about 1 millisecond and has 0.3 grams HMX base 
charge. To avoid any electrical noise and current the detonator wires 
were shielded and grounded. In addition the aluminum cylindrical shell 
was grounded, which helps in obtaining .. more clear oscilloscope traces. 

Several designs of plexiglass plugs were used and the final one used is 
shown in Fig. 22. The use of a plexiglass as a transmitter of the wave 

to the aluminum shell was advantageous in trying to obtain a symmetrical 
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Table 1. - Specimen Specifications 
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wave, since the wave in plexiglass travels slower than in the aluminum. 

. • / • * . 

For each explosion a new plug was used, since damage did occur to the plug 
after the explosion. 

A second type of electrical detonator used was 'the Du Pont X-8l 1 D 
miniature conductive mix detonator. This type of detonator has a faster 
rise time of 4-8 microseconds, and has 52 mg of explosive. There are no 
wires in this type of detonator and the problem of interferences was reduced. 
The fast firing time will enable us to obtain better oscilloscope traces. 
Since this type of detonator has more explosive, it was embedded in a lexan 
plug which is stronger and more resistant to impact. The basic design of 
the lexan plug was similar to the plexiglass plug. However, since the 
X-811D detonator are also smaller than the first type the lexan could be 
shortened and we could thus obtain a better impact condition at the 
aluminum shell. The lexan plug is shown in Fig. 23. 

Micromeasurement strain gages, 1/8 inch long, are mounted on the outside 
surface of the shell at various axial locations, as shown in Fig. 24. Ellis 
BAM-1B bridge amplifiers tuned to a frequency response of 100 KHZ are used 
to provide proper calibrated strain scaling of the oscilloscope. A trigger 
gage located near the impacted end was used to trigger the scope. Three 
alignment gages mounted on the outside surface 120 deg. apart were used to 
assure asimultaneity of impact, however in our testing only two gages were 
used. The arrival time of the propagating pulse at the various gage 
locations is obtained from the oscilloscope trace; since the distance between 
strain gages is known the velocity of propagation of the strain pulse is 
obtained by 
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The repeatability of the results was difficult since there appeared 
to be slight deviations between the detonators. However, the resulting 
strain traces have the same characteristics and are very similar in shape. 

A typical block diagram showing the circuitry, instruments, strain 
gages and the test specimens are shown in Fig. 25. 

RESULTS 

Typical oscilloscope traces are shown in Figs. 26 and 27. t Observation 
of these figures show the presence of a bending wave propagating at the 
plate velocity, Cp, which is faster than the shear velocity, c $ . Figure 26 
presents typical traces using the X-549 electrical detonator inserted in 
plexiglass. We obtained better results in this group of tests by using 
the circumferential gages no. 4 and 6 in Fig. 24. The measurement of the 
wave speeds yielded an average shear velocity of 

c s = 122,950 in/sec. 

The alignment gages show a reasonably good impact, but the presence of 
bending is shown by the difference in magnitude of the strain pulse between 
the alignment gages. However, a separation between the longitudinal and 
the shear wave are not clearly observed. Figure 27 presents typical traces 
using the X-811D electrical detonator embedded in lexan. Since the firing 
time is faster we used gages no. 7 and 9. The measurement of the wave 
speeds yielded 

c D = 214,500 in/sec. 

r 

c = 125,000 in/sec. 

These values are slightly higher than the analytical values of 

Cp = 205,333 . in/sec. 

c s = 113,338 in/ sec. 
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FIGUREZ5- TYPICAL BLOCK DIAGRAM OF CIRCUITRY , INSTRUMENTS 
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Observation of Fig. 27 s. hows that the results obtained from -the use of 
X-811D detonators are more consistent. The alignment gages show a better 
simultaneity of strain pulse than the previous types of electrical detonator. 

A better comparison between theoretical and experimental shear wave 
speeds is probably impossible, since, the analytically predicted value 
involves the use of a shear correction factor. The exact specification 
of this factor is questionable. 

CONCLUSION 

In this section an experimental technique for generating and measuring 
shear wave velocities in cylindrical shells under radial impact is presented. 
Although the strain pulses are not simple in shape the arrival time of the 
pulse at a gage location is shown very clearly, from which the wave speed, 
can be determined. There is still further investigation to be made in order 
to understand the transient response of a cylindrical shell to radial impact. 
However, the first step in this study was completed, namely the experimental 
verification of the speed of a shear wave propagating in cylindrical shell. 



VII. CONCLUSIONS AND RECOMMENDATIONS 


The conclusions of this study are briefly listed below. 

1. For clean cylindrical shells subjected to longitudinal impact: 

a. ) the single most important parameter effecting the magnitude. 

of the shell's transient response is the ratio h/x (x being 
related to the pulse duration, t q , by X = x Q Cp) . Except for 
the rectangular pulse, the peak values of strain or stress 
will always depend upon this 'ratio, 

b. ) single-pulse dispersion (wave spreading) is dependent on h/x 

In much the same manner as in the dispersion of continuous 
harmonic dispersion, 

c. ) the shape of the pulse behind the wave front becomes indist- 

inguishable for long pulse durations (small values of h/x) and 
all the shapes studied here become similar as the pulse travels 
down the shell , 

d. l except for the uniaxial theory, which predicts constant shape 

and magnitude of propagating pulse, the membrane (with and 
without shear) and bending theories predict, for practical 
purposes, similar transient responses, 

e. } the use of a sine series to approximate the pulse shapes is 

a useful tool for predicting conceptual or actual transient 
responses. The response for each of the sine terms is 
computed and the result for a particular shape is obtained 
by superposition of these individual results. 

2, For cylindrical shells having geometrical discontinuities subjected 
to longitudinal impact: 

a.) for short pulses having a short rise time (almost a step) 
the simple uniaxial result of 
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can be used to predict the effect of the geometrical .■ 
discontinuity on response. •- < 

b.) For longer pulse durations, the peak values occur behind 
the initial wave front and a computer calculation incor- 
porating matching conditions at the geometrical discontin- 
uity is required to predict the effect of the discontinuity 
on the transient response. 

3. The generation of a shear wave in a cylindrical shell has been 
accomplished with the measured shear velocity agreeing, within 8%, 'with the 
analytically predicted value. 

Recommendations for further work include. 

1. Experimental verification of the variation of strain response. with 
the h/x ratio, 

2. More analytical studies involving a larger range of geometrical 
discontinuities in the cylindrical shell, 

3. Experimental verification of the results for the shell having a 
geometrical discontinuity, 

4. An analytical study into the relationship between single pulse dis- 
persion and the classical harmonic dispersion, 

5. Further experimental work in the area of shear wave propagation and 
the transient response of shells to shear loadings, 

6. An analytical study into the effect of the shape of a shear pulse 
on a she! 1 's response, 

7. A joint analytical - experimental study into the response of a 
structure to a combined longitudinal - shear excitation. 
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VIII. GENERAL 


VIII. 1 - Relevance of Results 

The results of the study presented in this Report can be applied directly 
to the following technical areas: 

• a. Response of Structures to Pyrotechnic Loads 

Utilizing the results of thi s study into effects of pulse 
shapes or transient responses, we are in a position to predict 
the important parameters of the resulting loading pulse caused 
by a pyrotechnic event. In addition, we can apply the results 
of this study to the prediction of the effect of a structure's 
geometrical discontinuities on the resulting response of a 
structure subjected to a pyrotechnic event, 
b. Understanding Ultrasonic Signals in Structures (NDT) 

Since the thrust of this grant's study was the use of sine 
"building blocks", we have used the results of the transient 
response of the shell structure to these sine pulses in the 
interpretation of transmitted ultrasonic signals. The, single most 
important result used for this purpose is the dispersive effect 
caused by the thickness to wavelength ratio (h/x). For example, 
when transmitting an ultrasonic signal along a shell the magnitude 
(and shape) of this signal changes due to geometrical or material 
defects and dispersion. To date, in the field of NDT only limited 
information is available which will help in understanding which 
portion of the change in signal shape is due to defect and which 
portion is due to dispersion. The results of this year's grant 
are being applied to an ultrasonic study which will hopefully 
lead to a better understanding of this signal interpretation. The 
initial results are encouraging. 
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c, Response of Composite Material Structures to Dynamic Loads 

Much of the information and technology developed in this 
year's grant {and previous years') is in the process of being 
applied to composite material structures, both homogeneous.' and 
laminated. The initial results show that these capabilities 
are directly applicable to problems, such as, FOD in engine or 

helicopter blades and the dynamic loading of composite structures 

\ 

by foreign objects, pyrotechnic events , or docking. For example, 
one of our initial results have shown that for certain laminate 
plate configurations subjected to dynamic moment inputs, large 
normal and shear stresses are developed in the laminate, which 
is not the case in isotropic plates. Realizing that laminates 
cannot withstand high shear stresses, this result is seen to be 
important for the proper design of laminate structures to with- 
stand dynamic loads. 

VIII. 2 - Personnel Involvement 
3 - Faculty; Drs. R. Mortimer, P. Chou, and J. Rose 
3 - Graduate Students; Messrs. Blum, Rodini, and Cokonis 
2 - Undergraduate Students; Messrs. FI is and Nga Le 
1 - Technician; Mr. R. Tschirschnitz 

Of the three graduate students listed above, two are Ph.D. students 
and one is a Master's student. Mr. Blum has recently completed his Ph.D. 
thesis, which was supported by this grant. 
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IX. ADDITIONAL WORK COMPLETED UNDER GRANT 

In addition to the studies described in Sections II through VIII, 
other work related to the thrust of this Grant, and supported by this 

Grant, was completed. A brief description of thesetasks follows: 

♦ 

a. Equations Governing the Axisymmetric Motions of a Laminated 
Composite Cylindrical Shell. \ 

A system of dynamic equations of motion which govern the 
motions of a laminated composite cylindrical shell have been 
derived. This system of equations was derived by combining the 
usual isotropic shell theory development (Ref. 1) with the 
technique developed for anisotropic laminated plates (Ref. 10). 

This theory includes the effects of transverse shear, rotary 
and radial inertia, and bending. The transverse shear effect 
is extremely important when considering laminated structures. 

The resulting system of equations can be used to analyze the 
dynamic (transient) axisymmetric response of a cylindrical shell 
which is constructed from a number of isotropic or anisotropic 
layers. The form of the resulting system of equations is 

u" - — + ( ) v” = 

y 

V" - — + ( ) u" = 

y 

* x " - — + ( ) * ■ ■ - (ix-i) 

X 2 9 

2 

- — + ( ) - - - - - 

9 r 2 x 

z 

m W 

W H — =----- 

c 2 
3 
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where u, v, w, >p v and ip. denote the longitudinal, circumferential, 
and radial displacement, and the longitudinal and circumferential 
rotations, respectively, c-| , C2» 03 are wave speed parameters. The 
parenthesis appearing on the left-hand side of these equations 
involve material and geometric properties; the right-hand side of 
these equations involve the five displacements and their first 
derivatives. This theory can be, extended to include conical shells. 
Computer Code for the Dynamic Analysis of Laminated Structures. 

A computer code, based on the method of characteristics, 
capable of solving systems of equations equivalent to eq. (IX-1) 
has been written. Example runs to test the accuracy of this program 
still remain to be completed. This code, in addition to existing 
codes at Drexel , enable us to analyze the transient behavior of 
composite laminate shells, plates, and beams where each individual 
lamina may be isotropic or anisotropic. 

Cylindrical Shell with Finite Length Geometrical Discontinuity. 

y ' 

A computer code capable of analyzing the transient behavior 
of a shell having a finite length geometrical discontinuity has 
been developed. This code is an extension of the work reported in 
Ref. 5. 

Shear Wave Computer Code. 

A computer code capable of analyzing the transient response 
of structures to step loadings in transverse shear has been 
developed. 
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X. APPENDICIES 


APPENDIX A 


EQUATIONS GOVERNING THE MOTION OF A CYLINDRICAL SHELL 


Four systems of governing differential equations were used in this 
study. These equations, as presented here, are nondimensi onalized with 
respect to h such that 


\ . t c 

u T _ . - _ w - _ x ' - - _E 


u = fi’’^ = ^’ w = F ,x = 


T - 


The first system of equations includes membrane, bending, transverse 
shear, and rotary inertia effects and is given by (see Ref. 1) 
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The third system of equations represents the classical membrane theory 
and is given by . ' 

3 2 u _ a 2 u h 3w ■< 

(A-3) 

\ 

The final system of equations is the simple uniaxial theory and is given 
by 
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APPENDIX B 


LEAST SQUARES LINEAR ESTIMATION 

< i . 

Capabilities of Program 
The program will: 

1. Compute coefficients to any linear relation 

-\ 

2. Print and plot results on printout, 
e.g. 

F(x,y,z)^ A Q + A-j" V-j (x,y,z) + A 2 V 2 (x s y,z) ...A n V p 

F(x,y,z) - Function being estimated 
V-j(x,y»z)“- Estimating functions 
(x,y,z) - Independent coordinate variables 
A -»'A n — Coefficients computed by program 

N ote : 

(A) Approximation takes place on a finite interval 

(B) All functions must be bounded 

(C) All functions must be defined on the same coordinate system - 
Coordinate system can be comprised of any number of independent 
variables, e.g.: 

F (x) , F (x, r, e), F (z, Y, V ] , V 2 ) 

CPI All functions can be discontinuous". ■■ ' 
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SUoROGl INS F I L ‘‘‘AT ( 5,oS ) 

..U0U3LE.. Pi SCI S.L;:.t;__E.^fcLB,.AS. 

.1 V 0 I.-4& ;v S I ON 9(90,20) ,30(20,20) 

.9.0 • ■ __C ONMCN/ A P/ MO ,..N V ,NCV ,.N0.I ,.M0 1L IC 

91 NT=?!D*KV 


. . 9 . 2 . .... . 

.. ... MT.l-hHO . .... . . 

*' 

‘ ■ } 

93 

0 4 

NTL=NT+NDV 
I F ( N 0 - 1 ) 1 , ! . 2 



9 5 

2 CONTINUE 



_ 9.6 

. 0 031.4 M= 1 , MQ ... 




9 7 NP=NV 

98 : _..DJQ3J.5 K= 2.,AQ 

99 00301 1=1, N V 


101 

102. 

IF(E( M,L ) ) 312, 313, 312 
312 B(H,NP)=BlfM)«K 



103 

109 

GO TG 301 
313 B ( M . NP ) =0 . 



105 

106 

301 CONTINUE 

315 CONTINUE • 



107 

108 

314 CONTINUE 
• 1 CONTINUE 



109 

_L10_ 

IFIICI11, 10,11 

11. CONTINUE 



111 

112 

88 ( 1, 1)=N0 
00304K=1,NTL 



113 
1 14 

K S=K + 1 

. " A S = 0 . 



115 

116 

00305 M= 1 , NO 
AS=AS+B(M.K) 



117 

LL8 

305 CONTINUE •' 

B fi ( 1 . K S ) = AS 



119 

120 

B8(KS,1)=AS 
304 CONTINUE 



121 

122 

10 CONTINUE 
D0302N= 1 • NT 



123 
_L2 4 

■ N 1 = N + 1 C 
D0302 K=N . NTL 



125 

126 

K S = K+ I C 
Ab=0 • 



127 

128 

00303 M=1 , NO 
AS=B(M.K)*B(M.N)+AS 



129 

130 

3-03 CONTINUE _ - - - - — 

1 r ( NT 1 — K ) 4 » 3 « 3 

- - - • • 


131 

132 

3 CONTINUE 
B R ( Ml » KS ) - AS 



1 33 
134 

B 3 ( K S , N 1 ) =AS 
GO TO 302 



135 

_n '3 

1 } 7 
i. o c 

4 CONTINUE 
b B ( M , ix S ) = A S 

I'.Z COOT I : U 
-\ E f u .< o 


- 

139 

END 





in. 


.} z 


231 . 
23 3 
? * c 


240 

2tl 

242 

.24.3. 

244 

'223- 

246 

247 


250 


...5.UBR.0UL I.NE _ 50L V.( A.). 1 

0 0 U 2L if P 2 2 C I S I ; . A , X f A S ?,SA P 


i. ! -’EM 3 TO 


- (2 7, 


102 


C ■ i' ! • J • v / - ** / f ; 0 * ' '■< V * f i L. 7 y N - i I t l\u I 1 ! 

±_I.C ; 

m A i = :•! A + 1 

J:',l A. =.H 

.4 A L = i'i A + 1 ; D V 

J 10,1.0.0. L = LA . 

L 1 = L+ 1 

NS = -L0. . . 

D010 1 iN= L ? MA 

J.EJLAiK*JL LL1 02^10 1 ,..1.0.2 : 

NS = N 

GO TO 103 


101 CONTINUE 
10 3 I F ( NS- L ) 113 . 106.105 
248 105 D0104N=1 , HAL 

223 MUljji). 


1 04 


A ( L i N ) = A ( NS ♦ N ) 


252 

253 

• 10 6 

ASP=1./A(L,L) 
D0200N=L 1 , MAL 



254 


A ( L , N ) = A S F* A ( L , N ) 



25 5.. 

20 0 

CONTINUE 



256 


0 01 07 N=L 1 , M A 



257 


IF ( A (N, L) ) 108.107,108 



258 

103 

SAP=1 ./A (N,L ) 



259 


DG109K=L 1 , MAL 



260 


A (N,K)=A(N,K)*SAP -A ( L , K ) 



261 

109 

CONTINUE 



262 

107 

CONTINUE 



26 3 

100 

CONTINUE 



264 


00 IK* MAI, MAL 



265 


AlMAyK ) =A ( M A » K .. ) /'A { M A , M A ) 



266 

1 

CONTINUE 



267 

113 

CONTINUE 



268 


RETURN 



26 9 


END 



— 


— — ' -- - -• - - •• - ... . . .. .. . •. 

.... , . 


270 


SUBROUTINE DIAG(A,MA) 



271 


DOUBLE PRECISION A 



272 


DIMENSION A (20, 20) 



273 


C OMMON/ AP/MOyNV , NDV , ND I , ND , IC 



274 

f 

MA 1-MA+ 1 



275 


Ml A=MA-1 



276 


MAL=MA+NDV 



277 


DO 1 1 L=MA1 , MAL 



278 


D01 10N= 1 , Ml A 



279 


N 1 =MA 1-M 



200 


N 2 = N 1 - 1 



40 1 

1 1 4 

DO 111 K.~ 1 , NR 





•• (<\,L j = AU. ,i. ) - -- ( N , r i i. :: \ ( Vi , L ! 



c. 3 3 

ill 

CONTINUE 



2o^ 

110 

CONTINUE 



2o5 

1 1 . 

CONTINUE 



23 6 


RETURN ‘ • 



2 6 7 


END 




Y 


.IS*. 
16 5 

1 L 

1-77 


1 o 9 
170 


5 Uc '7 0 U T I £ P P I N T S ( A , B , C , M A ) 

OOUbLc- PRECISION A 4 B»P M » T Y , RR 

£ I -fe *1 S I • £. (20 ,20 ) , b (-0 , 20.) ,CUOOO ) 

C J 1**1 r**. ■)••!/ -P/ s'. 1 5, > 4 V ,NOV, ho I ’» NC » I C 
COMMON/ PP/VAR (20 , 10 ) » V - D ( 20 ,10 ) 


. C 0 i*i M ON A N G ( 9 ) , 6 L A N K , S T A P. 

ICO FORMAT! 1H 1 , / / » 2 OX * * E S T I 14 A T 

11,77) ' ' “7" 


P-. fi 


UK C T I • 0 . N 


= * , 20 7 . 


171 

44 

FORMA T< ///, 15X, * ESTIMATED • ,-10X , * CO NB I NED FUNCTIONS * ,9X , • 0 I 

FFER ANC * 


1 

, 1 4X , * X * ) 


172 

45 

FORMAT (// ,20X , • LINEARLY COMBINED ESTIMATING FUNCTIONS*,/ 

) 

173 

43 

FORMAT ( 4 ( 1 0 X , F 1 4 . 7 ) .) 


174 

40 

FORMAT ( 2CX , * ( * , L 14 .7 , 1 )*( * » 20A1, ')**(*, I 2, * )' ) 


175 

41 

F ORMAT ( / ' RMS DIFF = *,014.7',/) 


176 

.42 

FORMAT (• *,D14.7, r = CONSTANT*) •„ 


177 

4 6 

FORMAT (//,• COMPARISON BETWEEN ESTIMATED FUNCTION AND LIN 

EARLY CO 


1 

MBINED ESTIMATING FUNCTIONS') 


178 


K S = 1 - I C - - 


1 7 9 


MA1=MA + 1 


180 


M 1 A=MA- 1 


-18-1 


MA| =MA + NDV 


182 


NT = NV* ND 


-183- 


nti=nt+i 


184 


KK = (MA-IO/NV 


1 q 5 . 

00501 = MA1 .MAI 1 ..... 


136 


LL=NT1+L-MA 1 


1 37 


K = 1 -MA 


188 


PRINT 100 , (VAD( I ,K) , 1=1 ,20 ) 


1 89 


PRINT 45 


100 


I F( IC-1 )4, 2,4 


1 0 1 

2 

PR J MT4 7 f A ( 1 , 1 I 


192 

4 

CONTINUE 


193 


i = ir. 


194 


DQ1 I 0=1 , KK 


195 


nmiv=i.Nv '' 


196 


1 = 1 + 1 


19 7_ 

1 

PRINT 40 . A ( T .L ) . (VAR ( A , IV ) . J=1 ,70 ) . 10 


.198 


PRINT 46 


19.9— 


-E.8.I1LT .44 _ .. .... 


200 


RM=0. 


20 1 


DO 310N-1.N0 


202 


IF(KS-1 )11, 10,11 


-2.0.3 

1 0 

CONTINUE * 


204 


TY = A(1,L ) * B ( fl , K S ) 


-20.5- 


on to 1 7 


206 

11 

CONTINUE 


70 7 


T Y= A ( 1 . L ) 


208 

12 

C QMT I NU E 


..2.09— 


nrm 1 k = i .mi a 


210 


K 1 = K + 1 


211 


K2=<+KS 


7 > 

C. i. c. 

“5 "* 1 

TY = Br!,;<2)"A(Kl,L ) +TY 

t r;*j£ 


214 


0 ’. = T Y- 5 { , L L ) 


.21 5. . 


i_?.=i\n+N 


215 


L 3 = N 0*+ L 2 






PRtMT 5 

?ife C ( L 2 ) =f: ( N t LL i ' • _ 

21? - .CXI 2. l.=IY — - - • - 

2 1 0 C ( L 4)=RR 

.2.2.1.. _P.PXMI_.43 ,.a.( : ! U-L..J. ,T_Y_,?Xl.X i.M J. 

222 RM=RR*RR+-F.ri 

Z23 3.LC.XGNXI.N U E : ' ; 

224 RM=RM«* . 5 /FLO AT {MO ) 

2 ? 5 WXIXEJjZrALlBii : 

22 6 K=L-MA 

227 CXLl^.P-LIiXlJ<i..C^riUjA..OJ_ 

228 50 CONTINUE 

2 29 RETURN 

230 END 


140 SUBROUTINE SORT! A » BB »MA> 


141 

142 

DOUBLE PRECISION A,B3,BS 
DIMENSION A (20,20) ,BB< 20,20) 


143 

COMMON/ AP/NDrNV*NOV»NDI f NDt IC 


144 

MA1=MA+1 


145 

M IA=MA- 1 ^ - 


146 

MAL=MA+NDV 


147 

NT=ND*NV+IC 


148 

NT1=NT+ 1 


.149 

00307 1 = 1 , Ml A 


150 

1 1=1 + 1 


151 

A ( I , I ) = 38 ( 1,1) 


152 

00307 11 , MA 


153 

BS=BB(I,J) 


154 

A ( I , J ) = 8 S 


155 

A ( J , I ) = 3 S 


156 

307 CONTINUE 


157 

" A ( MA » M A } = B B { M A , MA ) 


158 

0011=1, HA 



159 D01L=MA1 r MAL 

16 0 L L = Nil + Lj-_M AL 

161 1 All ,L>=bB< I. ,LL J 

162 RET URN 


END 




.C 

SU’j.-: .UTiiic iPLCT 



r 

c 

PUkPlSH - 



c 

c 

PLOT SEVERAL CRGSS-V AM I-ABL ES VcRSUS A oASL V AA I A Hi £ 
DIRECT PLOT I NO OF INPUT VARIABLES (ML- INTER PCLAT I CM ) 


c 

c 

USAGE 



c 

c 

CALL APLOT(NO,A,N,M,NS) 



c 

c 

DESCRIPTION OF PARAMETERS 

NO - CHART NUMBER (3 DIGITS MAXIMUM) 



c 

c 

A - MATRIX OF DATA TO BE PLOTTED. FIRST .COLUMN REPRESENTS 
BASE VARIABLE AND SUCCESSIVE COLUMNS' ARE THE CROSS- 


c 

c 

c 

o. 

VARIABLES (MAXIMUM IS 9). 

N - NUMBER OF ROWS IN MATRIX A 
M - NUMBER OF COLUMNS IN MATRIX A (EQUAL TO 
MIIMP.FR "HF VAR! ARIFS). MAXIMUM IS. 10.. 

THE TOTAL 


c 

r 

NS -CODE FOR SORTING THE BASE VARIABLE DATA 
ORDER 

IN ASCENDING 


c 

r 

0 SORTING IS NOT NECESSARY (ALREADY 
ORDER). 

IN ASCENDING 


C 

r. 

1 SORTING IS NECESSARY. 



c 

r. 

REMARKS 

THF FOl LOW f MG FXFf IITT ON CARDS NEEDED TM MAIN PROGRAM 



c 

r 

COMMON ANG (9 ), BLANK, STAR 

RFAD I 8 . SI Rl A NK, (A NG (I). T=1.9). STAR 



C 

c 

5 FORMAT (11A1) 

THF FOLLOWING DATA CARD NFFDFD IN MAIN PROGRAM 

i 


c 

c... 

STARTING IN COULUMN #1 

.173456789* 



c 

c 

4 ’ IS A BLANK SPACE 



c 

c 

SUBROUTINES AMD FUNCTION SUBPROGRAMS REQUIRED 
NOME 



c 

c 


t 

\ 


c 


4 

.2.3 A. 


SUBROUTINE API DT ( NO , A , N, M, NS ) 


2 89 
2 9Q 


DIMENSION OUT (116) , Y P R ( 1 1 ) , A(1000) 

x- COMMON ANG ( 9 ) ,RI ANK, STAR ~ — “T‘V.Z: 

. - 

791 

c 

1 F0RMATDH1.6nx.7H CHART f I3,//) 


292 

293.. 


2F0RMAT(1H »F11.4,5X,1Q1A1 ) 

3 FORMAT ( / . POX. 1 >>>>>>>>>> ESTIMATING F 

u it C T I 0 N s 

2° ■“ 

. . ..... 

1 <<<<<<<<<<*,//) 

7 FORMAT (in , 1 6 X , 1 C 1 H . 

> 



L • * • * * 

F JEVTTd- '.ov fl icir. 4) 



c 

. ...c 

c 

V " 



.L9 tu. 

...... ;;ll=.i . . . . . . .. . 



* 

5 ? r - 

‘ 



22X 


APLOT C cont 'j 

L 


C 



299 00 14 J — i , i-. 

■30.0 . I.F..LA U. ) - A ( J ) 1 • 1 A 1—14 » 1.1 

501 11 .1=1 —N 


303 


— 

LL-J'lt __ - .. 

DO 12 K = 1,M 


— 


30 4 



L = 1 +N 




305 



LL=LL+N 




306 



F = A ( L ) - 




307 



A ( L ) = A ( L L ) 




308 


12 

A ( L L ) = F 

•• 



309 


14 

CONTINUE 




3.10. . 


_L5 

CONTINUE 





C 







C 


TEST NLL 





c . 






311 


_I6 

IF(NLL) 20, 18, 2 C 




312 

. c 

18 

NLL- 50 





o o 


PRINT TITLE 




313 


20 

WRITE (6, 1)N0 




314 



IF(NO.EQ.0)PRINT3 





c 







c 


FIND SCALF FOR CROSS-VARIABLES 





G 
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IF(A( JJ-YMAX) 40,40,30 
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Y M I N= A ( J ) 
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GO TO 40 
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YMAX= A( J ) 
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CONTINUE 
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